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Abstract 

A critical point symmetry for the prolate to oblate shape phase transition is introduced, 
starting from the Bohr Hamiltonian and approximately separating variables for 7 = 30°. 
Parameter- free (up to overall scale factors) predictions for spectra and B(E2) transition 
rates are found to be in good agreement with experimental data for 194 Pt, which is supposed 
to be located very close to the prolate to oblate critical point, as well as for its neighbours 
( 192 Pt, 196 Pt). 
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1. Introduction 

Critical point symmetries in nuclear structure are recently receiving considerable at- 
tention [1, 2, 3], since they provide parameter-free (up to overall scale factors) predictions 
supported by experimental evidence [4, 5, 6, 7]. So far the E(5) [U(5) (vibrational) to 0(6) 
(7-unstable)] [1, 4, 5] and the X(5) [U(5) to SU(3) (prolate deformed)] [2, 6, 7] critical 
point symmetries have been considered, with the recent addition of Y(5) [3], related to 
the transition from axial to triaxial shapes. All these critical point symmetries have been 
constructed by considering the original Bohr equation [8], separating the collective (3 and 
7 variables, and making different assumpions about the u(/3) and 7/(7) potentials involved. 

Furthermore, it has been demonstrated [9] that experimental data in the Hf-Hg mass 
region indicate the presence of a prolate to oblate shape phase transition, the nucleus 194 Pt 
being the closest one to the critical point. No critical point symmetry for the prolate to 
oblate shape phase transition originating from the Bohr equation has been given so far, 
although it has been suggested [10, 11] that the (parameter-dependent) 0(6) limit of the 
Interacting Boson Model (IBM) [12] can serve as the critical point of this transition, since 
various physical quantities exhibit a drastic change of behaviour at 0(6), as they should 
[13]. 

In the present work a parameter-free (up to overall scale factors) critical point symmetry, 
to be called Z(5), is introduced for the prolate to oblate shape phase transition, leading to 
parameter-free predictions which compare very well with the experimental data for 194 Pt. 
The path followed for constructing the Z(5) critical point symmetry is described here: 

1) Separation of variables in the Bohr equation [8] is achieved by assuming 7 = 30°. 
When considering the transition from 7 = 0° (prolate) to 7 = 60° (oblate), it is reasonable 
to expect that the triaxial region (0° < 7 < 60°) will be crossed, 7 = 30° lying in its middle. 
Indeed, there is experimental evidence supporting this assumption [14]. 

2) For 7 = 30° the K quantum number (angular momentum projection on the body- 
fixed i'-axis) is not a good quantum number any more, but a, the angular momentum 
projection on the body-fixed x'-axis is, as found [15] in the study of the triaxial rotator 
[16, 17]. 

3) Assuming an infinite well potential in the /3-variable and a harmonic oscillator po- 
tential having a minimum at 7 = 30° in the 7- variable, the Z(5) model is obtained. 

On these choices, the following comments apply: 

1) Taking 7 = 30° does not mean that rigid triaxial shapes are prefered. In fact, it has 
been pointed out [18] that a nucleus in a 7-flat potential [19] (as it should be expected for a 
prolate to oblate shape phase transition) oscillates uniformly over 7 from 7 = 0° to 7 = 60°, 
having an average value of 7 a „ = 30°, and, therefore, the triaxial case to which it should 
be compared is the one with 7 = 30°. Furthermore, it is known [20] that many predictions 
of models involving large rigid triaxiality are very close to the predictions of 7-soft models 
involving 7-fluctuations such that •jrigid of the former equals r y rms of the latter. In addition, 
the equivalence between 7-instability and rigid triaxiality with 7 = 30° has been shown in 
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relation to the 0(6) limit of IBM using projection techniques [21, 22]. In view of these, it 
is not surprising that the Z(5) predictions describe well nuclei like 194 Pt, which are known 
to be good examples of the 0(6) symmetry [20]. 

2) Taking an infinite well potential in (3 (while 7 is fixed at 30°) corresponds to a 
transition from a triaxial vibrator to a triaxial rotator [15, 16, 17], in the same way that an 
infinite well potential in (3 in the X(5) model (in which 7 = 0° is assumed [2]) corresponds 
to a transition from a vibrator [U(5)] to a prolate rotator [SU(3)]. This point will be further 
discussed in Section 7. 

3) In view of the above, it is not surprising that the Z(5) model gives results compatible 
with earlier work on the prolate to oblate shape phase transition [9, 10]. In the earlier 
work [9, 10], the 7-soft pool in which the critical point of the prolate to oblate transition is 
expected to lie, is crossed by moving from 7 = 0° to 7 = 60° in the (3 7^ region, i.e. away 
from the vibrational ((3 = 0) regime. When moving from SU(3) (prolate) to SU(3) (oblate) 
on the appropriate side of the extended [10] Casten triangle [20], one then identifies 0(6) 
as the critical point. In the Z(5) model the same pool is crossed in a different way, by fixing 
7 = 30° and moving from the triaxial vibrator (close to f3 — 0) to the triaxial rotator (far 
from/3 = 0) [15, 16, 17]. 

In Sections 2 and 3 of the present work the /3-part and the 7-part of the spectrum will be 
considered respectively, while B(E2) transition rates will be studied in Section 4. Numerical 
results will be reported in Section 5 and compared to experimental data in Section 6, while 
Section 7 contains a summary of the present results and plans for further work. 

2. The /5-part of the spectrum 

The original Bohr Hamiltonian [8] is 



2B 



d 4 _g_ , 1 g • S — — V Q 

V f)R + R2 c ;„ Q A , tU, S111 A/12 2^ 



/3 4 d/T dd ' /3 2 sin3 7 97^"'97 4/3 2 ^ ^ ( 7 _ |^) +V ^^ 

(1) 

where (3 and 7 are the usual collective coordinates, while Qk (k = 1, 2, 3) are the components 
of angular momentum and B is the mass parameter. 

In the case in which the potential has a minimum around 7 = 7r/6 one can write the 
last term of Eq. (1) in the form 

£ 2( Ql 2~kT\ w $ + + Ql) = ml + QI + Ql) - 3Q?. (2) 

fc=lj2)3 sm 2 (7 - fk) 

Using this result in the Schrodinger equation corresponding to the Hamiltonian of Eq. (1), 
introducing [2] reduced energies e = 2BE/K 2 and reduced potentials u = 2BV/H 2 , and 
assuming [2] that the reduced potential can be separated into two terms, one depending on 
(3 and the other depending on 7, i.e. 7) = u((3) + ^(7), the Schrodinger equation can 
be separated into two equations 



ZL,M = ^L,a(/3), (3) 
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1 d d | ' 

(P 2 ) sin 37 c?7 



77(7) = ^77(7), 



(4) 



where L is the angular momentum quantum number, a is the projection of the angular 
momentum on the body-fixed £'-axis (a has to be an even integer [15]), (P 2 ) is the average 
of (3 2 over £(/?), and e — ep + e 7 . 

The total wave function should have the form 



*(/3,7,^)=a,a(/3)^(7)^M )Q (^), 



(5) 



where 0j (i = 1, 2, 3) are the Euler angles, denote Wigner functions of them, L are the 

eigenvalues of angular momentum, while M and a are the eigenvalues of the projections of 
angular momentum on the laboratory fixed i-axis and the body-fixed x'-axis respectively. 

Instead of the projection a of the angular momentum on the x'-axis, it is customary to 
introduce the wobbling quantum number [15, 23] n w = L — a. Inserting a = L — n w in Eq. 
(3) one obtains 



tL, nw {P)=tt)ZL,n w (P), (6) 



where the wobbling quantum number n w labels a series of bands with L = n w ,n w + 2,n w + 
4, . . . (with n w > 0) next to the ground state band (with n w = 0) [15]. 
In the case in which u((3) is an infinite well potential 



u(P) 



if P < Pw 
00 for P > p w 



(7) 



one can use the transformation [2] = /5 3//2 ^(/5), as well as the definitions [2] ep = kp, 
z = Pkp, in order to bring Eq. (6) into the form of a Bessel equation 



d 2 i 1 di 

— - H + 

dz 2 z dz 



with 



4L(L + 1) - 3a 2 + 9 JL(L + 4) + 3n w (2L - n 



2 2 
Then the boundary condition £(Pw) — determines the spectrum 

X Q 



t f3;s,v — £l3;s,n w ,L ~ (&s,i/) i h 



s,v 



s,v 



Pw ' 



and the eigenfunctions 

Z.AP) = ts,n w ,M = 6,a,L(/5) = C s , v p-*l 2 J v (k a , v P\ 



(8) 



(9) 



(10) 
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where x Sjl/ is the sth zero of the Bessel function J v (z), while the normalization constants 
c S}U are determined from the normalization condition / °° (3 4 ^ 2 v (f3)df3 = 1. The notation for 
the roots has been kept the same as in Ref. [2], while for the energies the notation E Sjnw>L 
will be used. The ground state band corresponds to s — 1, n w — 0. We shall refer to 
the model corresponding to this solution as Z(5) (which is not meant as a group label), in 
analogy to the E(5) [1], X(5) [2], and Y(5) [3] models. 
3. The 7-part of the spectrum 

The 7-part of the spectrum is obtained from Eq. (4), which can be simply rewritten as 



d 2 cos 37 d \ 



{(3 2 ) \d^ 2 sin 37 



77(7) = ^77(7). 



;i2) 



As already mentioned, we consider a harmonic oscillator potential having a minimum at 
7 = 7r/6, i.e. 

7T\ 2 



"(7) 



1 



\ c 7 



6/ 



1-2 ~ 71 

2 C7 ' 7 = 7 "6- 



(13) 



In the case of 7 ~ 7r/6 the cos 37 term vanishes and the above equation can be brought 
into the form 

d 2 1 



77(7) = 6^(^)77(7) 



#7 2 ' 2 

which is a simple harmonic oscillator equation with energy eigenvalues 



2c 



and eigenfunctions 

^(7) = iV n ,^(67)e^ 2/2 , 
with normalization constant 



0,1,2,... 



(14) 



(15) 



\ v /i2 n ^n^!' 

Similar potentials and solutions in the 7-variable have been considered in [8, 24] 
The total energy in the case of the Z(5) model is then 

E(s, n w , L, n^) = E + A(x SyV ) 2 + Bn^. 



(16) 



(17) 



(18) 



It should be noticed that in Eq. (14) there is a latent dependence on s, L, and n w 
"hidden" in the {(3 2 ) t erm. The approximate separation of the (5 and 7 variables is achieved 
by considering an adiabatic limit, as in the X(5) case [2, 25]. 

4. B(E2) transition rates 
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The quadrupole operator is given by 



cos (7 - |) + ^(PgW + 2>2_ 2 (*)) - (7 - |) 



(19) 



where t is a scale factor, while in the Wigner functions the quantum number a appears next 
to and the quantity 7 — 2n/3 in the trigonometric functions is obtained from 7 — 2nk/3 
for k — 1, since in the present case the projection a along the body- fixed x'-axis is used. 
For 7 ~ 7r/6 this expression is simplified into 

If 2) = --^^(^S(^) + 2>2-2(«i))- (20) 
B(E2) transition rates are given by 

where the reduced matrix element is obtained through the Wigner-Eckart theorem 

(L f a f \T^\L t a t ) = W^aA \T^\ \L m ). (22) 

V 2L / + 1 

The symmetrized wave function reads 



HP,lA) = ts, a ,Mrin, (7) \l 16 X(l + ) a ) + (-!) L ^U (23) 



where the normalization factor occurs from the standard integrals involving two Wigner 
functions [26] and is the same as in [15]. a has to be an even integer [15], while for a = 
it is clear that only even values of L are allowed, since the symmetrized wave function is 
vanishing otherwise. 

In the calculation of the matrix elements of Eq. (22) the integral over 7 leads to unity 
[because of the normalization of 77(7)], the integral over (3 takes the form 

((3)(3 4 df3, (24) 

where the (3 factor comes from Eq. (20), and the f3 A factor comes from the volume element 
[8], while the integral over the angles is calculated using the standard integrals involving 
three Wigner functions [26]. The separation of the integrals occurs because 77(7) does not 
depend on a, while in £(/3) only even values of a appear. The final result reads 

5 t 2 1 

B(E2; LiCti — > £/«/) 



16tt2 (l + 5 ai:0 )(l + 5 af , ) 
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(Li2L f \ai2a f ) + {LflL^ai - 2a f ) + (-l) Li (LfiL f \ - afiaf)] 2 7j( Si , L h on, s f , Lf, a f ). 

(25) 

One can easily see that the Clebsch-Gordan coefficients (CGCs) appearing in this equation 
impose a Aa = ±2 selection rule. Indeed, the first CGC is nonvanishing only if ai + 2 = af, 
while the second CGC is nonvanishing only if cej — 2 = af. The third CGC is nonvanishing 
only if ai + aj = 2, which can be valid only in a few special cases. The angular part of this 
equation is equivalent to the results obtained in [15]. 

The ground state band (gsb) is characterized by n w = L — a = 0. Therefore transitions 
within the gsb are characterized by ai = Li and af = Lf. Normalizing the B(E2) rates to 
the lowest transition within the gsb we obtain 

o „ , 9 n B(E2; (L + 2) g — > L g ) 5 2L + 1 7j(l, L + 2, L + 2, 1, L, L) 

R g , g {L + 2^L)- B{m _ ^ - ^ - g ^5(1 + d Lfl ) IJYJ^W) ' 

(26) 

where the (1 + 5l,o) factor comes from the fact that to all transitions within the gsb only the 
second CGC in Eq. (25) contributes, except for the lowest one, to which both the second 
and the third terms contribute. 

The even levels of the 7!-band are characterized by n w = L — a = 2, which means 
a = L — 2. Using the same normalization as above one obtains for the transitions from the 
even levels of the 71 band to the gsb 

o (T n B(E2;L,-,L 9 ) 15 Zj(l, L, L - 2, 1, L, L) 

7e " e " M ' B(E2;2 g ^0 g ) (L + 1)(2L + 3) 7j(l, 2, 2, 1, 0, 0) ' 

(27) 

where the (I + Sl^) factor is due to the fact that for all transitions only the first CGC of Eq. 
(25) contributes, except in the case of 2 7 — > 2 g , in which both the first and the third terms 
contribute. The angular parts of Eqs. (26) and (27) coincide with the results obtained in 
[27]. 

In a similar manner the following ratios are also derived 

p (T T 1 u B(E2;L 1 ^(L + l) g ) 5 I 2 (l, L, L - 1, 1, L + 1, L + 1) 
7odd ' sl + } B(E2;2 g ^0 g ) L + 2 P(l, 2, 2, 1, 0, 0) 1 j 



B(£2; 2, - fl ) 



5(2L-1)L(2L + 1) , 7 2 (l,L + 2,L,l,L,L-2) 

1 + M ^7- - n n , , (29) 



2(2L + 3)(L + 2)(2L + 5) V 7 2 (1, 2, 2, 1, 0, 0) 

R iodd^~todd( L + 2 — ► 7,) 



7?(£2; (L + 2) 7 - L 7 ) 



B(£2; 2 9 - a ) 

5L(2L + 1) 7 2 (1,L + 2,L + 1,1,L,L-1) 
2(L + 2)(2L + 5) 7 2 (1,2,2,1,0,0) ' 



(30) 
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5(2L-3)(2L- 1) 
L(L + 1)(2L + 1) 




(L-L-l) 



(1 + M 



B(E2-L 1 — > (L — 1) 7 ) 
5(£2;2 S ^0 9 ) 

I 2 (1,L,L-1,1,L- l,L-3) 
/ 2 (1,2,2,1,0,0) 



(31) 



It should be noticed that quadrupole moments vanish, because of the Aa = ±2 selection 
rule, since in the relevant matrix elements of the quadrupole operator one should have 

OLi = OLj. 

5. Numerical results 

The lowest bands of the Z(5) model are given in Table 1. The notation L SjTlw is used. All 
levels are measured from the ground state, O^o, and are normalized to the first excited state, 
2^0- The ground state band is characterized by s — 1, n w — 0, while the even and the odd 
levels of the 71-band are characterized by s — 1, n w — 2, and s — 1, n w — 1 respectively. 
The /Si-band is characterized by s = 2, n w = 0. All these bands are characterized by 
n^y = 0, and, as seen from Eq. (18), are parameter free. The fact that the 71-band is 
characterized by = is not surprising, since this is in general the case in the framework 
of the rot at ion- vibration model [28]. 

B(E2) transition rates, normalized to the one between the two lowest states, B(E2;2i )0 — > 
Oi^), are given in Table 2. 

6. Comparison to experiment 

Several energy levels and B(E2) transition rates predicted by the Z(5) model are com- 
pared in Table 3 to the corresponding experimental quantities of 194 Pt [29], which has been 
suggested [9] to lie very close to the prolate to oblate critical point. Its neighbours, 192 Pt 
[30] and 196 Pt [31], which demonstrate quite similar behaviour, are also shown. Not only 
the levels of the ground state band are well reproduced (below the backbending) , but in 
addition the bandheads of the 71-band and the /9i-band are very well reproduced, without 
involving any free parameter. The staggering of the theoretical levels within the 71-band 
is quite stronger than the one seen experimentally, as it is expected [18] for models related 
to the triaxial rotator [15, 16, 17]. 

The main features of the B(E2) transition rates are also well reproduced. As far as 
the transitions from the 71-band to the ground state band are concerned, the transitions 
£1,2 — ► £1,0 are strong, while the transitions (L + 2)i j2 — > £1,0, which are forbidden in the 
Z(5) framework, are weaker by two or three orders of magnitude. Even the augmentation 
of B(E2;2i )2 — > 2 lj0 ) relative to B(E2;4i j2 — > 4 lj0 ), which is due to a mathematical detail, as 
explained below Eq. (27), is very well seen experimentally. 

7. Discussion 

In summary, a critical point symmetry for the prolate to oblate shape phase transition 
has been introduced by approximately separating variables in the Bohr Hamiltonian for 
7 = 30°. The parameter free (up to overall scale factors) predictions of the model, called 
Z(5), are in good agreement with experimental data for 194 Pt, which is supposed to lie close 
to the prolate to oblate critical point [9], as well as for its neighbours ( 192 Pt, 196 Pt). 
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In addition to the points made in the introduction, the following comments apply: 

1) The /3-equation [Eq. (6)] obtained after approximately separating variables in the 
Bohr Hamiltonian is also exactly soluble [32, 33] when plugging in it the Davidson potentials 
[34] 

=P 2 + f 2 , (32) 

where (3q is the minimum of the potential. In analogy to earlier work in the E(5) and X(5) 
frameworks [35], it is expected that (3$ = should correspond to a triaxial vibrator, while 
Po — > oo should lead to a triaxial rotator [15, 16, 17]. 

2) Using the variational procedure developed recently in the E(5) and X(5) frameworks 
[35], one should be able to prove that the Z(5) model can be obtained from the Davidson 
potentials by maximizing the rate of change of various measures of collectivity with respect 
to the parameter (3 , thus proving that Z(5) is also the critical point symmetry of the 
transition from a triaxial vibrator to a triaxial rotator. 

Work in these directions is in progress. 
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Table 1: Energy levels of the Z(5) model (with = 0), measured from the L SjUw = O^o 
ground state and normalized to the 2 10 lowest excited state. See Section 5 for further 
details. 
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Table 2: B(E2) transition rates of the Z(5) model, normalized to the transition between 
the two lowest states, B(E2;2 lj0 0i j0 ). See Sections 4 and 5 for further details. 
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i4 l,0 


iz l,0 


o.^oo 


J-^1,2 


1 9, „ 
-1^1,2 


9 ROD 
z.uyu 


loi,i 


loi,i 


3 1 on 
0. iyu 


16i,o 


14i, 


3.642 


16l,2 


14l,2 


2.938 


17i,i 


15i,i 


3.392 


18i >0 


16i, 


3.803 


18i, 2 


16i )2 


3.151 


19i,i 


17i,i 


3.566 


20i i0 


18i, 


3.944 


20i i2 


18i, 2 


3.335 








2l,2 


2i,o 


1.620 


3i,i 


4i,o 


1.243 


3i,i 


2l,2 


2.171 


4l,2 


4i,o 


0.348 


5i,i 


61,0 


0.972 


5i,i 


4l,2 


1.313 


6l,2 


61,0 


0.198 


7i,i 


81,0 


0.808 


7i,i 


6l,2 


1.260 


8l,2 


81,0 


0.129 


9i,i 


10i, 


0.696 


9i,i 


8l,2 


1.164 


10l )2 


10i, 


0.092 


Hi,i 


12i, 


0.614 


Hi,i 


10l, 2 


1.069 


12l )2 


12i, 


0.069 


13i,i 


14i, 


0.551 


13i,i 


12i, 2 


0.984 


14l,2 


14i, 


0.054 


15i,i 


16i, 


0.507 


15i,i 


14l,2 


0.910 


16l )2 


16i, 


0.043 


17i,i 


18i, 


0.459 


17i,i 


16l, 2 


0.846 


18l,2 


18i, 


0.035 


19i,i 


20i, 


0.425 


19i,i 


18i, 2 


0.790 


20i, 2 


20i )0 


0.030 
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Table 3: Comparison of the Z(5) predictions for energy levels (left part) and B(E2) tran- 
sition rates (right part) to experimental data for 192 Pt [30], 194 Pt [29], and 196 Pt [31]. See 
Section 6 for further discussion. 





Z(5) 


192p t 


194p t 


196p t 


L ii) 


L (f) 


Z(5) 


192p t 


194p t 


196p t 


4i,o 


2.350 


2.479 


2.470 


2.465 


4i,o 


2i,o 


1.590 


1.559 


1.724 


1.476 


61,0 


3.984 


4.314 


4.299 


4.290 










81,0 


5.877 


6.377 


6.392 


6.333 


4l,2 


2l,2 


0.736 




0.446 


0.715 


10i )0 


8.019 


8.624 




8.558 


6l,2 


4l,2 


1.031 






1.208 


2l,2 


1.837 


1.935 


1.894 


1.936 


3i,i 


2l,2 


2.171 


1.786 






4l,2 


4.420 


3.795 


3.743 


3.636 














6l,2 


7.063 


5.905 


5.863 


5.644 


2l,2 


0i, 


0.000 


0.009 


0.006 


0.0004 












2l,2 


2i,o 


1.620 


1.909 


1.805 




3i,i 


2.597 


2.910 


2.809 


2.854 


4l,2 


2i,o 


0.000 




0.004 


0.014 


5i,i 


4.634 


4.682 


4.563 


4.526 


4l,2 


4i,o 


0.348 




0.406 




7i,i 


6.869 


6.677 






6l,2 


4i,o 


0.000 






0.012 


2 ,o 


3.913 


3.776 


3.858 


3.944 
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